Accepted for Publication in The Astrophysical Journal 

Preprint typeset using I4TgX style emulateapj v. 10/09/06 



o 
o 

(N 
D 



o 
u 

6 



FULL POLARIZATION SPECTRA OF 3C 279 
D. C. Homan', M. L. Lister\ H. D. Aller\ M. F. Aller\ J. F. C. Wardle'* 

Accepted for Publication in The Astrophysical Journal 

ABSTRACT 

We report the results of parsec-scale, multi-frequency VLBA observations of the core region of 3C 279 in 
Stokes /, linear polarization, and circular polarization. These full polarization spectra are modeled by radiative 
transfer simulations to constrain the magnetic field and particle properties of the parsec-scale jet in 3C279. 
We find that the polarization properties of the core region, including the amount of linear polarization, the 
amount and sign of Faraday rotation, and the amount and sign of circular polarization can be explained by a 
consistent physical picture. The base of the jet, component D, is modeled as an inhomogeneous Blandford- 
Konigl style conical jet dominated by a vector-ordered poloidal magnetic field along the jet axis, and we 
estimate its net magnetic flux. This poloidal field is responsible for the linear and circular polarization from this 
inhomogeneous component. Farther down the jet the magnetic field in two homogeneous features is dominated 
by local shocks and a smaller fraction of vector-ordered poloidal field remains along the jet axis. This remaining 
poloidal field provides internal Faraday rotation which drives Faraday conversion of linear polarization into 
circular polarization from these components. In this picture, we find the jet to be kinetically dominated by 
protons with the radiating particles being dominated by electrons at an approximate fraction of > 75%, still 
allowing the potential for a significant admixture of positrons. Based on the amounts of Faraday conversion 
deduced for the homogeneous components, we find a plausible range for the lower cutoff in the relativistic 
particle energy spectrum to be 5 < 7; < 35. The physical picture described here is not unique if the observed 
Faraday rotation and depolarization occur in screens external to the jet; however, we find the joint explanation 
of linear and circular polarization observations from a single set of magnetic fields and particle properties 
internal to the jet to be compelling evidence for this picture. 

Subject headings: galaxies : active — galaxies: jets — galaxies: individual: 3C 279 — radiation mechanisms: 
non-thermal — radio continuum: galaxies 
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1. INTRODUCTION 

The three-dimensional magnetic field structures and the 
particle populations of extragalactic jets from Active Galactic 
Nuclei (AG N) are still not we ll understood. iLaing & Bridld 
(l2002h and iLaing et aP (l2006h have recently made progress 
in studying the 3-D magnetic field structures of kilo-parsec 
scale jets, but little is known about the 3-D field structure 
near the jet origin, on parsec or sub-parsec scales. We 
wish to know if the jet magnetic field shows a structure 
which has its roots in the magnetic field in the supermassive 
black-hole/accretion-disk system responsible for giving rise 
to the iets (e.g.[B landford & Znaiek 1977; Meier et al. 2001- 
iKoide et alJl2002l;IVlahakis & Koni gl 2004). iMarscher et all 
(I2002h has suggested a direct disk-jet connection on the ba- 
sis of X-ray/radio correlations, and this connection may ex- 
tend to the magnetic field which threads the accretion disk 
around the black hole. Unanswered questions include, for 
example, is there significant vector-ordered poloidal field 
along the jet axis or perhaps a toroidal/helical field struc- 
ture confining the jet and indicating a je t current jAsada et alj 
I2002t iGabuzda et alJl2004l: lFendtll2006l) ? It is also unknown 
whether the particle population of jets is prima rily electron- 
proton, p^e~, or electron-positron, e'^e~ (e.g. IWardle et al] 
[l998;.Sikora & Madeiski 2000; De Youngi2006h . We param- 
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eterize this unknown as the lepton number, £ = (ri- - n+) / (n_ + 
n+), where «_ and «+ are the number densities of electrons 
and positrons respectively. Additionally, the limits on the 
power-law particle spectrum which gives rise to the observed 
synchrotron radiation are poorly constrained. The relativistic 
number density can be parameterized as N^d'^ = K^'^d^ for 
7/ < 7 < 7» where it is assumed to have a hard cutoff to the 
power-law at both high and low energies. ICelotti & FabianI 
(1993) found that the lower cutoff 7/ set the scale for the 
bulk kinetic luminosity of jets because the low energy par- 
ticles dominate the particle density. 

With sub-milliarcsecond resolution. Very Long Baseline 
Array (VLBA) observations of extragalactic radio jets can 
study synchrotron emission from jets within a few parsecs 
of the jet origin. Measurements of linear polarization from 
milli-arcsecond jets are sensitive to the net magnetic field or- 
der in the plane of the sky; however, multiple 3-D magnetic 
field structures can all yield the same observed linear polar- 
ization. For example, a transverse shock of tangled mag- 
netic field, a toroidal field, or a high-pitched helical field will 
all produce linear polarization with the electric field vector 
parallel to the jet axis. Likewise, a vector-ordered poloidal 
field, shear in a tangled field, or a low-pitch helical field will 
all produce linear polarization with the electric field vector 
transverse to the jet axis. Faraday rotation of linear polar- 
ization is sensitive to magnetic field order along the line of 
sight and to the properties of the particles doing the rotation. 
The observed Faraday rotation must, in some cases, result 
from thermal particles and magnetic fields which are along 
the line of sight but external to the jet, but for rotations smaller 
than about 45°, it is difficult to distinguish internal fro m ex- 
ternal Faraday rotation without additional information jBurnl 
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Il966l) . Recently, several authors have attributed Faraday ro- 
tation gradients observed approximately transverse to jets as 
evidence of toroidal or helical magnetic fields either in the jet 
itself or in a s heath l ayer around the j et (Asada et al. 2002; 
Qabuzda et al. 2004; Zavala & Tayloi] [20051: lAttridge et alJ 
2005; Gomez et al. 2008); however, these gradients are of- 
ten gradients in magnitude of a single sign of the observed 
rotation and do not show the clear anti-symmetric signa- 
ture expected for such fields, so an additional external Fara- 
day screen must be proposed as well. Alternatively the ob- 
served gradients in rotation measure could be due to den- 
sity gradients in the material surrounding the jet. In the case 
of 3C 273, the case for toroidal or helical fields is stronger 
as the same dire ction of the gradient is observed at multi- 
ple jet loca tions (Zavala & Tavlor 2005; Attridge et al. 2005|; 
I Asada etal.. 2008). Marscher et al. (2008) have also argued 
for helical field structure in BL Lacertae based on temporal 
rotations in linear polarization angle from the jet core region. 

Parsec-scale observations of circular polarization from 
extra-galactic jets provide additional constraints that can 
break some of the degeneracies inherent in linear polariza- 
tion observations. Circular polarization may be generated ei- 
ther as an intrinsic component of the emitted synchrotron ra- 
diation or via Faraday co nversion of Unear polarization into 
circular polarization (e.g. Jones 1988). Like Faraday rotation, 
conversion is a bi-refringence effect; however, unlike rotation, 
conversion is much stronger in relativistic particles than ther- 
mal particles, and we don't expect significant Faraday con- 
version from magnetic fields and particles external to the jet 
(iJones & O'Dellll 19771: iHoman et al.ll2001h . In this way, cir- 
cular polarization probes the jet magnetic fields and particles 
directly without modification from external screens. 

Parsec-scale circular polarization observations of AGN jets 
have b e en reported b y Wardle et al. ( 1998); Homan & Wardle 
dl 999^; 'Homan et al. (2001); Homan & Wardle (2003, 2004); 
[Homan & Lister (2006), and most recently by Gabuzda et al. 
12008). Most AGN jets appear to have less than 0. 1 to 0.2% 
circular polarization, with 10% to 20% of jets detected at the 
level of ~ 0.3% to 1 .0% of the Stokes / emission from or very 
near the base of the jet (or "jet core"). The highest levels of 
circular polarization detected are 2 to 4% of the local Stokes 
/ emission in the co re region of the nearby radio galaxy 3C 84 
(iHoman & Wardl e|l2b04) and in the in tra-day variable source 
PKS 1519-273 (Macqua rt et al.ll2000l) . Single dish monitor- 
ing of circular polarization from AGN jets by the UMRAO has 
been ongoing since 2002 ( Aller et al. 2003), and the Austrian 
Compact Telescope Arra y has studied the in tegrated circular 
polarization from AGN (Rayner et al.ll2000h . Circular polar- 
ization has also been observed in intra-day variable sources, 
micro-quasar s, low luminosity AGN, and the Galactic Center 
(iBower et al.l[T99 9. 2002; Brunthaler et al. 2001; Fender et al. 
120001 120021; iMacquart et al.ll2000t iSault & M acquart 19991) . 
To date, no strong correlations have been found between the 
appea rance of circular polarization and other source proper- 
ties (iHoman et alJl200lt iRayner et al.ll2000t iHom an & ListeJ 
120061) 7 however, the lack of strong correlations may simply 
be related to the small fraction of detected sources and the 
low levels of circular polarization in those objects. There 
is increasing evidence that at least some circularly polarized 
sources tend to have a preferred "handedness" of circular 
polarization (jKomesa roff et al. 1984; Homan & Wardle 1999; 
IHoman et a P 2001: lower et al . 2002; Homan & Lister 2006) 
suggesting a persistent magnetic field structure is responsible 
for setting the sign of the observed polarization, although it is 



important to not e that changes in sign have been observed in 
some cases re.g. iAlleretal.ll2003l) . 

To date, radiative transfer modeling of parsec-scale cir- 
cular polarization has been at only one or two frequencies, 
limiting our abihty to uni quely constrain physical jet prop- 
erties. IWardle et al] (IT998) preferred a model for 3C279 
where circular polarization was produced via Faraday con- 
version of linear into circular polarization and this conver- 
sion was driven by some internal Faraday rotation. Their re- 
sults implied a low cutoff in the relativistic particle spectrum 
of 7/ < 20, suggesting the jet was dominated by electron- 
positron pairs to keep the jet's kinetic lumin osity within 
reasonable bounds (e. g., IC elotti & Fab ianlll993h . How ever, 
subsequent work by [Rus zkowski & B egeknanl (l2002h and 
iBeckert & Falckd (1200^ showed that some addition of ther- 
mal matter to the jet allowed a larger value for 7/, reducing the 
need f or a e lectron-positron dominate d jet. iHoman & Wardld 
( 120031) and iHoman & Wardld ( |2004 also preferred Faraday 
conversion mo d els for PKS 0607-157 and 3C 84 respectively. 
iGabuzda et al.l (HoOSj) have linked the observation of trans- 
verse rotation measure gradients with parsec-scale circular 
polarization observations in a qualitative fashion to show that 
both observations are consistent with helical magnetic fields if 
all eight jets in their sample emerge from the south magnetic 
pole of the central engine. 

To make progress on using parsec-scale circular and lin- 
ear polarization observations of AGN jets to constrain their 
magnetic field structure and particle populations, we have em- 
barked on a program to study the linear polarization, circular 
polarization, and Stokes / spectra of AGN jets at six frequen- 
cies with the VLBA. We call these "full polarization spectra" 
of AGN jets, as we use all four Stokes parameters to constrain 
radiative transfer models. We combine these single-epoch, 
parsec-scale data with integrated monitoring by the UMRAO 
at 4.8, 8.0, and 14.5 GHz. The first VLBA observations in this 
program took place in November of 2005 for 18 jets. We de- 
tected strong circular polarization at multiple frequencies in 
three jets: 3C 84, 3C 279, and 3C 380. Here we report our re- 
sults on 3C 279. Future papers will explore 3C 84 and 3C 380, 
as well as the sources in our sample with little or no detected 
circular polarization. In ^we describe our observations, cal- 
ibration, and Gaussian modeling of the core region of 3C 279 
to obtain spectra for fitting. In ^we describe a variety of pos- 
sible models for producing circular and linear polarization in 
jets and then use analytic and computational radiative trans- 
fer models to fit the full polarization spectra obtained in ^ 
We discuss the results of the radiative transfer modeling in 
H and our conclusions appear in ^ Throughout the paper 
we assume a cosmology where Hg = 70 km/s/Mpc, fi^ = 0.3, 
and ilA = 0.7, and so at the redshift of 3C279 (z = 0.538), one 
milli-arcsecond corresponds to a projected linear scale of 6.34 
parsecs. 

2. OBSERVATIONS 

2.1. VLBA Observations and Reduction 

In 2005, November 17-20, we observed 18 AGN radio jets 
with the National Radio Astronomy Observatory's Very Long 
Baseline Array (VLBA)^ at six frequencies: 8.01 GHz, 8.81 
GHz, 12.35 GHz, 15.37 GHz, 22.23 GHz, and 24.35 GHz. 
The observations were recorded at each frequency in dual cir- 

The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
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cular polarization with single-bit recording at four intermedi- 
ate frequencies (IFs) of 8 MHz bandwidth for a total band- 
width of 32 MHz. The observations were made over a contin- 
uous 72 hour period to minimize the effects of source variabil- 
ity with each of three 24 hour segments devoted to a pair of 
frequencies: 8.0 and 22 GHz, 12 and 15 GHz, and 8.8 and 24 
GHz respectively. Each source was visited 9 times at each ob- 
serving frequency with these "scans" highly interleaved with 
neighboring sources to maximize (u,v)-plane and parallactic 
angle coverage. Scan lengths were approximately 300 sec- 
onds at 22 and 24 GHz, 215 seconds at 12 and 15 GHz, and 
130 seconds at 8.0 and 8.8 GHz. Each source thus obtained 
approximately 45 minutes integration time per frequency at 
22 and 24 GHz, 32 minutes at 12 and 15 GHz, and 19 minutes 
at 8.0 and 8.8 GHz. Our calibrated Stokes I and polarization 
images at each of these frequency bands appear in Figures 1, 
2, and 3 respectively. 

A-priori data calibration, fringe-fitting, self-calibration, and 
calibration for linear and circular polarization were performed 
at Deni son University using the techniques described in detail 
by [Liste r & Homan (2005) and Homan & Lister (2006). At 
1 5 .4 GHz, the electric vector position angle (EVPA) of the lin- 
ear polarization was calibrated by comparin g to the MOJAVE 
database of rotated D-terms, as described in iLister & HomanI 
(|2005). At the other frequencies, we used the polarization of 
eight well defined jet components and looked at their devia- 
tions from the 15 GHz calibration (which was taken as ab- 
solute); although there was some clear Faraday Rotation in 
three of the components, the offsets at each frequency were 
remarkably consistent. We did make corrections for the ap- 
parent rotation measures of those three components, but our 
calibration is not strongly dependent on those values. Indeed, 
if we rely on only the five components without apparent Fara- 
day rotation, we get the same calibration to within a degree 
at each frequency. Based on the scatter of the offsets between 
components, we estimate that our EVPA calibration is good 
to better than one degree at 12 GHz and below, and good to 
better than two degrees at 22 and 24 GHz. These numbers 
are relative to the 15 GHz calib ration, which we expec t to be 
accurate to within one degree, (iLister & Homanll2005l) : how- 
ever, relative calibration between the frequencies is our main 
concern here for the modeling of the emission region. In our 
results we therefore consider only calibration uncertainty rela- 
tive to 15 GHz added in quadrature with modeling uncertainty. 

Circular polarization calibration was d one via the gain 
transf er technique described in detail by iHoman & Listed 
(|2006), including their Monte Carlo methods to estimate the 
final uncertainty in the observed circular polarization. We 
have made some small improvements upon the methods de- 
scribed there to allow semi-automatic flagging of antennas 
and some individual scans that show a high RCP/LCP gain 
ratio. In particular we flagged an antenna if the standard devi- 
ation of its RCP/LCP antenna gain ratio was more than twice 
the median standard deviation of all the antennas. This re- 
quirement only resulted in the flagging of the OV antenna at 
8.0 GHz which had a very large (factor of ~ 2) gain discrep- 
ancy between the right and left hand feeds. We also flagged 
individual scans if their RCP/LCP ratio deviated by more than 
5 times the standard deviation from the smoothed gains. In 
practice this accounted for a very small percentage, 0.1 -0.3 
percent, of the overal l scans being flagged. As described in 
iHoman & Lister! (1200 6). we also excluded sources with low 
SNR, as indicated by large RCP/LCP gain fluctuations, or 
with apparently high levels of circular polarization, as indi- 



TABLE 1 

Circular Polarization of Core Region. 



Measured Values 


If All Stokes-V on "D" 


Freq. 


V 




(GHz) 


(mJy) 


(%) 


(1) 


(2) 


(3) 


8.01 


50.8 ±9.6 


1.99 ±0.38 


8.81 


48.2±8.0 


1.55 ±0.26 


12.35 






15.37 


47.6±7.9 


0.88 ±0.15 


22.23 


58.2± 13.4 


0.86 ±0.20 


24.35 


71.0± 16.5 


0.94 ±0.22 


Note. 


— Column,s are 


as follows: (1) Frequency of 



observation in GHz; (2) Stokes-V flux density for the core 
region in mJy; (3) Fractional circular polarization under 
tile liypotlietical scenario, discussed in §3.4.1, where all 
the Stokes-V flux is associated with coinponent "D" in the 
core region. 

cated by a larger RCP/LCP gain offset from the smoothed 
gains, from contributing to the final smoothed gain-transfer 
table which would determine the overall circular polarization 
calibration.* These requirements excluded only between one 
and three sources from contributing to the final gain smooth- 
ing at each frequency, except at 24 GHz where seven sources 
were excluded, likely due to lower overall SNR at that fre- 
quency. Our circular polarization observations at 12 GHz 
were corrupted by strong radio frequency interference (RFI), 
likely from geostationary satellites passing within the tele- 
scope beam, so we do not report circular polarization results 
at that frequency. 

As discussed in detail in lHoman & List er (2006), phase cal- 
ibration for circular polarization is complicated if a source 
has significant extended structure in Stokes /. This may lead 
to spurious anti-symmetric structure in the resulting circularly 
polarized images due to phase errors in earlier steps in the cal- 
ibration. To account for this effect, we use the same approach 
described and tested by lHoman & Listed (12006) of adding an 
extra round of phase self-calibration, assuming no circular po- 
larization in the data. In general this extra round of phase 
self-calibration is very effective in removing genuine phase 
gain errors whil e preserving the o riginal circularly polarized 
signal; although iHoman & Listed (l2006i) did find that the am- 
plitude of the circularly polarized signal may be reduced by 
a few percent up to 10% and the position of the circular po- 
larization was shifted a small amount, less than about half a 
beam-width, to better align with the source peak. Real cir- 
cularly polarized structure significantly away from the source 
peak was not shifted by th e procedure. In the case of a source 
with a broad core region, iHoman & Listed (l2006h found that 
the circular polarization may be spread by this procedure to 
encompass the whole region. 

In recognition of the potential uncertainty of the precise 
placement of the circular polarization within the core region 
of 3C279, we measured the integrated amount of circular po- 
larization in the core region at each frequency using a single 
Gaussian component in the (u,v)-plane, and these results are 
reported in Table 1 . 

As described below, the core region of 3C279 is well mod- 

* The formal limit on the RCP/LCP gain fluctuations for low SNR sources 
was if the source had a standard deviation larger than twice the median stan- 
dard deviation of all sources. The formal limit for a source with large ap- 
parent circular polarization was 0.5%, as indicated by a systematic RCP/LCP 
gain offset of > 0.005 froin the smoothed gains. 
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(b) 8.8 GHz 
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Relative R.A. (marcsec) 



Fig. 1 . — Contour images of Stokes-I flux density (left), linearly polarized flux density (middle) and circularly polarized flux density (right) at 8.0 GHz (panel 
(a)) and 8.8 GHz (panel (b)). Contour levels begin at 5 mJy/beam for Stokes-I and circular polarization, and 10 mjy/beam for linear polarization, and increase 
in steps of x V2. Tick marks on the linear polarization image represent the measured electric vector position angle. A single contour from the Stokes-I image 
bounds the linear and circular polarization images to show registration. A cross-figure representing the FWHM dimensions of the restoring beam appears in the 
lower left hand comer of each panel. 
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(a) 12.3 GHz 




Relative R.A. (marcsec) 



(b) 15.4 GHz 




^5 -10 -15 -20 

Relative R.A. (marcsec) 

Fig. 2. — Contour images of Stokes-I flux density (left), linearly polarized flux density (middle) and circularly polarized flux density (right) at 12.3 GHz 
(panel (a)) and 15.4 GHz (panel (b)). Contour levels begin at 5 mJy/beam for Stokes-I and circular polarization, and 10 mJy/beam for lineal' polarization, and 
increase in steps of X \/2. Tick marks on the Hnear polarization image represent the measured electric vector position angle. A single contour from the Stokes-I 
image bounds the linear and circular polarization images to show registration. Note that severe RFI at 12.3 GHz prevented reliable measurement of the circular 
polarization at this frequency, so those data are not presented, see 0for a more detailed description of this problem. A cross-figure representing the FWHM 
dimensions of the restoring beam appears in the lower left hand corner of each panel. 
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Fig. 3. — Contour images of Stoke.s-I flux density (left), linearly polarized flux density (middle) and circularly polarized flux density (right) at 22.2 GHz 
(panel (a)) and 24.4 GHz (panel (b)). Contour levels begin at 10 mJy/beam and increase in steps of X \/7.. Tick marks on the linear polarization image represent 
the measured electric vector position angle. A single contour from the Stokes-I image bounds the linear and circular polarization images to shovi' registration. A 
cross-figure representing the FWHM dimensions of the restoring beam appears in the lower left hand comer of each panel. 
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1 -1 

Relative R.A. (marcsec) 

Fig. 4. — Uniformly weighted contour image of tlie core region at 24 GHz 
witli the main core components labeled. Note that component positions and 
FWHM sizes plotted here are from a single possible set of model components. 
As described in the text, we fit seven possible sets of components to span 
a range of plausible scenarios for fitting the structure of this source. For 
example, the two weaker, unlabeled components in this figure are fit as a 
single component in two of our seven possible models. 

eled in Stokes / and linear polarization by three closely spaced 
components labeled D, 5, and 4 in Figure 4, and the third 
column of Table 1 shows the fractional circular polarization 
at each frequency if all of the measured circular polarization 
originated from component D. However, as described above, 
the complications of phase calibration for circular polariza- 
tion limit our ability to confidently divide the measured cir- 
cular polarization between these components at the lower fre- 
quencies, and indeed at 8.0 and 8.8 GHz, the measured cir- 
cular polarization could come from any combination of these 
three components. At 15 GHz, the circular polarization is 
most likely associated with some combination of D and 5. At 
22 and 24 GHz, we are confidently able to assign the mea- 
sured circular polarization to component D at the base of the 
jet. To confirm this, we ran tests of our phase calibration pro- 
cedure at 22 and 24 GHz by generating simulated data with 
the same Stokes / structure as the clean component models at 
these frequencies, but with the addition of a 10 or 20 mJy cir- 
cularly polarized component at the locations of components 5 
or 4 with the remaining circular polarization on component D. 
We found that the extra round of phase-calibration assuming 
zero circular polarization was not able to move circular polar- 
ization that originated on components 5 or 4 to component D. 
These results give us confidence that the circular polarization 
observed on component D indeed originated on component D 
and was not transferred there by our phase calibration. 

2.2. Gaussian Model of the Core Region 



TABLE 2 

Flux and Linear Polarization of Core Region 
Components. 



Component 


Freq. 


/ 


mi 


X 


ID 


(GHz) 


(Jy) 


(%) 


(degrees) 


(1) 


(2) 


(3) 


(4) 


(5) 




8.01 


7 -i- 101 
z jjj ZLL yj.jyjL 


3 3 ± 7 


1 18 9 ± 4 




8 81 




3 9 ± 5 


^ j3 J -j- J 7 




12.35 


A ftf^O-l- 441 


3 Q_{_ Q 3 


1 18 7 ± 2 7 




15.37 




2.5 ± 0.2 


1 12.6 ±3.1 




22.23 


U. / jy ^1 \J.'~rJ 1 


2 4-1- Q 1 


5 ± 3 1 




24.35 




2.2 ±0.1 


1 (If) Q -L 2 s 


5 


8.01 


3.914± 0.355 


5.4±0.8 


18.1±5.4 




8.81 


3.865 ±0.354 


7.5 ±0.9 


19.9± 1.9 




12.35 


3.944 ±0.331 


6.8 ±0.8 


37.4 ±2.4 




15.37 


3.329 ±0.246 


8.6±0.6 


41.6± 1.8 




22.23 


2.812 ±0.236 


10.0 ±0.9 


49.3 ±2.2 




24.35 


2.977 ±0.260 


9.0±0.8 


48.8 ±2.0 


4 


8.01 


2.964 ±0.433 


13.3±1.7 


163.1 ±3.3 




8.81 


2.817 ±0.397 


13.5±1.3 


164.4 ±2.1 




12.35 


3.178±0.313 


13.5±0.6 


174.3 ± 1.5 




15.37 


2.668 ±0.228 


16.3±1.0 


175.6± 1.2 




22.23 


2.235 ±0.157 


16.6±0.8 


181.2±2.2 




24.35 


2.211±0.160 


16.3±0.8 


179.3 ± 1.9 



Note. — Columns are as follows: (I) Component identifier; (2) Frequency 
of observation in GHz; (3) Stokes-I flux density in Jy; (4) Percent fractional 
linear polarization; (5) Electric vector position angle of linear polarization. 

We used the Caltech VLBI program Difmap (IShepherdI 
Il997h to model the Stokes / structure and linear polarization 
of 3C279 at each frequency with particular attention to the 
structure of the core region which produces the observed cir- 
cular polarization. One challenge to this process is that the 
angular resolution of our observations differs by a factor of 
three between 8.0 and 24 GHz. To fit a comparable model 
over this range of resolution, we required that the relative sep- 
aration, sizes, and shapes of the Gaussian components remain 
the same from one frequency to another; however, we allowed 
the fluxes of each component to vary as well as the overall 
phase center of the model. 

One difficulty with this approach is that no single collec- 
tion of model components will be ideal for all frequencies, so 
in total we repeated this process to fit seven separate sets of 
model components in this fashion. Each of these separate sets 
of model components explored somewhat different assump- 
tions for how the structure of the core region and the remain- 
der of the source was fit. For example, in some cases we fit 
elliptical Gaussians to the core region and in other cases cir- 
cular Gaussian. We also tried varying the relative positions 
of the core components as well as the properties and num- 
bers of components used to describe the rest of the jet. Taken 
together, these seven different sets of model components sam- 
pled a range of reasonable possibilities for fitting the structure 
of 3C279. Overall the results for the core region were quite 
similar from one possibility to the next, and in Table 2 we 
report the average fluxes and polarizations from these differ- 
ent possibilities for each of the main core components. The 
uncertainties for the component properties given in Table 2 in- 
clude calibration uncertainties^ added in quadrature with the 
standard deviation of the component values from the seven 
different sets of model components. Figure 4 shows the size 
and spacing of the three core region components from a rep- 
resentative set of model components. 

^ The overall flux calibration at each frequency was assumed to be go od to 
5% based on the results reported in the appendix of lHoman et alj i2002l) . 
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Fig. 5. — Intergrated polarization monitoring by the UMRAO at 4.8, 8.0. 
and 14.5 GHz from 2003 through 2007. From the bottom to top. Stokes I, 
fractional linear polarization, electric vector position angle, and fractional 
circular polarization. Two week averages of the data are shown. The stan- 
dard errors, dominated by random measurement errors, are often smaller than 
the plotting symbols used. The horizontal line marks a zero level for circular 
polarization, and the vertical line indicates the epoch of the VLBA measure- 
ments. 

2.3. UMRAO Integrated Monitoring of3C279 

The integrated total flux density, linear polarization, and cir- 
cular polarization of 3C 279 have been monitored with the 
University of Michigan 26 meter telescope operating alter- 
nately at 4.8, 8.0 and 14.5 GHz. The prime-focus polarime- 
ters, utilizing rotating quarter-wave plates feeding orthogo- 
nal linearly-polarized transducers, measure all four Stokes 
parameters simultaneously. The general observing and cali - 
bration procedures used are described in lAUer et aTl (l2003h . 
Each series of on-off polarization observations, lasting ap- 
proximately 40 minutes, was preceded by position scans to 
verify the pointing of the telescope. Observations of 3C 
279 were interleaved with observations of reference sources 
at roughly two-hour intervals. The instrumental polarization 
was checked from observations of bright galactic HII regions 
which are assumed to be unpolarized. All observations are 
restricted to within three hours of the meridian to minimize 
possible instrumental effects at large hour angles. Sources 
are not observed within 15 degrees of the sun at 14.5 or 8.0 
GHz or within 30 degrees of the sun at 4.8 GHz to avoid solar 
interference. Because of these restrictions, there are annual 
gaps in the data for 3C 279 during September through mid- 
November at 4.8 GHz and from late September through late 
October at the two higher frequencies. Figure |5] shows the 
integrated emission of 3C 279 from 2003 through 2007. The 
annual gaps in the data are due to the close proximity of the 
sun at those times. We detect non-zero circular polarization at 
all three observing frequencies, and within the measurement 
uncertainties the Michigan instrument measured the same am- 
plitude and polarity of circular polarization as observed by the 
VLBA. During the period shown, circular polarization exhib- 
ited a preference to be negative at 4.8 GHz and positive at the 
two higher frequencies. 



3. DERIVING PHYSICAL PARAMETERS OF THE EMISSION 
REGIONS 

Our objective in this section is to model the magnetic field 
structure and particle properties in the core region of the 
parsec-scale jet of 3C279. Tables 1 and 2 provide a large 
number of constraints: Stokes / spectra for each component, 
fractional linear polarization (LP), nii, as a function of fre- 
quency for each component, electric vector position angle 
(EVPA), X, as a function of frequency for each component, 
and finally, fractional circular polarization (CP), nic, for the 
entire core region at low frequency and for component "D" at 
high frequency. 

In the analysis that follows, we take the jet to have a bulk 
Lorentz factor of F = 15 at an angle of 9 = 1.5° to the line 
of sight, giving a Doppler factor of 5 = 26 for the jet. These 
n umbers are consisten t with the recent kinematical analysis 
of iHoman et al.l (l2003h . Simil ar but slightly larger values are 
found bv 'J orstad et all (|2004|) : F = 20 and 6* = 1.0°, giving 
S = 36. Both sets of values give an observed angle in the 
frame of the emitting fluid of 9' ~ 40°, and our results are 
not strongly dependent on which set of values we choose. If 
we did choose the "larger" value for the Doppler factor, it be- 
comes somewhat easier to produce circular polarization via 
the intrinsic mechanism, as the estimated field strength s de - 
rived for the emission regions are somewhat larger (see 33.3b . 

3.1. Major Models of CP Production 

Before we begin the detailed analysis, we will briefly out- 
line five major models of circular polari zation production 
which may be relevant to this source (e.g. lWardle & HomanI 
I2003L and refer ences therein). 

Model (1) is stochastic production of circular polarization 
from a purely tangled magnetic field. In this model, the CP 
could be produced by either the intrinsic or Faraday con- 
version mechanisms, or some combination of the two. This 
model requires a fairly coarse-grained magnetic field tangling 
to do the job, as the sign of the produced circular polarization 
will differ from cell-to-cell, leading to square-root-N style 
cancellation. A key prediction of this model is that the circular 
polarization should vary in sign and amount across different 
frequencies, which probe different parts of the jet near opti- 
cal depth equals unity. We would also expect in this model 
that the CP would vary with time as the details of the mag- 
netic field tangling change as the jet flows outward. Our ob- 
servations of 3C 279 clearly show a consistency of sign and 
magnitude of circular polarization across the frequency bands, 
and previous work (e.g. lHoman & Listeri2006i) has shown that 
3C 279 has maintained the same sign of circular polarization 
at high frequency for a decade now despite having a fast flow- 
ing jet and several outbursts over that period. The UMRAO 
has also observed a consistent sign of positive CP in 3C 279 
at 14.5 GHz from 2003 through 2007, see figure |5l however, 
the UMRAO also finds periodic switches to negative CP at 
4.8 GHz which we discuss in ^ as a possible opacity effect. 
Based on these results, we conclude that stochastic production 
of CP in a random magnetic field is not a major mechanism 
producing the CP we see in 3C279. In the analysis that fol- 
lows, we will therefore assume that any tangled portion of the 
magnetic field in 3C279 is tangled on a very short length- 
scale to minimize the production of stochastic CP and allow 
exploration of the contribution by organized magnetic fields. 

Model (2) is intrinsic circular polarization from a strong, 
vector-ordered magnetic field, hereafter intrinsic circular po- 
larization. Given the consistency across epoch and frequency 
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in the observed circular polarization of 3C 279, the most plau- 
sible candidate for such a vector-ordered field in the jet is a 
poloidal field along the jet axis. We can estimate the strength 
of such a field. If the magnetic field were completely uni- 
form (no reversals or tangling), the expected circular polar- 
ization would be of order m^, ~ \/vb/ ^emit where vb = 2.8B 
MHz is the electron gyro-freq uency and Vemit = ^^^(1 + z)/S 
is the emitted frequency (e.g., IWardle & Homanll2003h . So 
for 1 % CP at 22 GHz, the required magnetic field strength is 
about B ^ 50 mG. This is not an unreasonable field strength, 
and as we shall see in §3.4.2, is of the order of the estimated 
field strength in component "D". 

It is important to note that the above estimate assumes not 
only a uniform magnetic field, but also that all of the emitting 
particles are electrons. This will be the case in a "normal" 
matter jet consisting of a pure electron-proton plasma. Allow- 
ing for the possibility of some admixture of electron-positron 
pairs, parameterized by the lepton number, £, defined in the 
introduction, the fractional circular polarization is nic cx £ and 
the required field strength would scale as B oc £~^. 

Model (3) is Faraday conversion of linear polarization into 
circular polarization driven by Faraday rotation, hereafter ro- 
tation driven conversion. In this model the LP is produced by 
whatever magnetic field order is available in the jet, such as 
shocked or sheared magnetic field. However, that LP cannot 
be directly converted into CP by the same field order because 
there needs to be some offset between the position angle of 
the LP and the angle of the magnetic field doing the conver- 
sion. The angular offset is provided by internal Faraday ro- 
tation within the jet, and in this model, the Faraday rotation 
depth is relatively small so that the field at the front of the jet 
is converting the LP emitted by the back of the jet into some 
CP. 

For this to be a consistent (not stochastic) process, the Fara- 
day rotation must be provided by some vector-ordered field 
in the jet, probably a poloidal field along the jet axis as in 
model (2). Additionally, there must also be some preponder- 
ance of electrons in the jet relative to positrons, i.e., I > Q. 
However, the requirements for some vector-ordered field and 
some preponderance of electrons in the jet plasma are not 
nearly as strong as they are for model (2), where the CP is 
produced entirely by the intrinsic mechanism. Model (3) has 
the additional requirement of low energy particles within the 
jet to produce the internal Faraday rotation, either due to a 
low cutoff, 7/, in the relativistic particle power-law spectrum 
N^d^ = K^~''d^ for 7/ < 7 < 7,, or due to the addition of some 
"cold", non-relativistic thermal matter to the jet. 

It is interesting to note that the similarity in requirements 
between models (2) and (3) means that they will often act 
together to at least some degree. 

Model (4) is a high rotation depth version of model 
(3), proposed and investigated in dependently by 
iRuszkowski & BegelmanI ([2002) and iBeckert & Falckd 
(I2002h . The requirements are almost the same as for model 
(3), but the rotation depth is much larger, so large, in fact, 
that significant rotation and conversion can happen over 
very small length scales in the jet. As long as this sense 
of rotation is common in all parts of the jet, no large scale 
magnetic field order is required (except for the vector-ordered 
field which produces the rotation) and net circular polar- 
ization can be produced in potenti ally large amount s with 
little if any net linear polariza tion dBeckert & Falckd 120021: 
IRuszkowski & Begelmanll2002l) . 



Model (5) is Faraday conversion from a helical magnetic 
field or other field structure which varies systematically in 
orientation across the jet. In this model, no internal Fara- 
day rotation is required, as the field orientation at the back 
of the jet is already at some angle with respect to the field 
orientation at the fro nt of the jet (e.g. lWardle & Homanll2003l : 
iGabuzda et aT]|2008l) . 

3.2. Numerical Modeling 

We have written a numerical simulat ion that solves the full- 
Stoke s equations of radiative transfer (^ Jones & 0'DelllfT977l: 
lJonesi ri988) either by numerical integratio n or, optionally, 
by usi ng the exact solutions calculated by i Jones & O'DeUl 
(Il977l). An early version of this simulation was used by 
IWardle et al.l (Il998h to interpret the first parsec-scale circu- 
lar polarization observations on 3C 279. The simulation mod- 
els the emission from a "homogeneous" line of sight broken 
into cells where each cell has the same physical properties, 
including spectral index, a, lepton number, £, Doppler fac- 
tor, 5, viewing angle, 6, and low energy power-law cutoff for 
the relativistic particle distribution, 7;. The exception to this 
is the magnetic field, which may consist of multiple compo- 
nents which can vary in direction and/or magnitude from cell 
to cell depending on the magnetic field model applied. Addi- 
tionally, multiple lines of sight can be put tog ether to construct 
simple inhomogeneous models (see ^3.4.2l i or homogeneous 
models where the magnetic field has transverse structure, such 
as toroidal or helical magnetic fields (see ^3.4. lb . By adjust- 
ing the physical parameters and magnetic field model applied, 
we can simulate any of the five conceptual models described 
in the previous section. The radiative transfer is done in the 
frame co-moving with the fluid, and here we have assumed 
a single flow velocity where the relativistic flow parameters 
are derived from observed pattern speeds for this jet given by 
VLBA kinematical analyses as described in the introduction 
to this section. This simulation is static in the sense that the 
magnetic field and particle properties of the emission region 
are assumed not to change during a light crossing time. 

3.2. L Parameterization of the Magnetic Field 

The magnetic field in our simulation can have three differ- 
ent components: a vector-ordered uniform field along the jet 
axis, Bi, = X /,„ a toroidal field, B, = Bj, x f, x (p/ pjet), and 
a randomly ordered field which varies stochastically from cell 
to cell, Br = B-i,x (l-fu-ft)- Here B* is a scaling factor which 
allows us to match a desired average perpendicular magnetic 
field strength, < B± >, set by observation, /„ and /, param- 
eterize the degree of uniform and toroidal fields respectively, 
and p/ Pjet is the fractional distance of a cell from the cen- 
ter axis of the jet, where pj Pjet = 1.0 is a cell at the outside 
edge of the jet. This parameterization of the toroidal field 
assumes that the current carried by the jet is uniformly dis- 
tributed. Note that non-zero values for /„ and /, together will 
produce a helical field order 

In addition to the three field components described above, 
the magnetic field may also be shocked, by shortening unit 
length to length, k< 1.0. This shock is assumed to be a trans- 
verse shock, so the shortening occurs along the jet axis. The 
result is amplification of the magnetic field components trans- 
verse to the jet axis as described in detail in the appendix of 
IWardle etali([T99l . 

3.3. Modeling the Stokes I Spectra 
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Fig. 6. — Stokes / spectrum of the core region of 3C 279. Component D 
is plotted witli red open squares. Components 5 and 4 are plotted witli blue 
X symbols and green filled triangles respectively. The combined core region 
flux is plotted in black filled squares. An estimate for the combined core 
region flux at 4.8 GHz from comparison of the UM RAO and VLBA data is 
plotted as an open star. The spectral fits described in j|3.3l are plotted as solid 
lines of the same color as the corresponding components. 

Figure |6] shows the Stokes / spectra of components D, 5, 
and 4 and their total flux density along with analytical models 
for each component. The total core flux at 4.8 GHz is esti- 
mated by comparing quasi-simultaneous UMRAO integrated 
measurements at 4.8, 8.0, and 14.5 GHz to our VLBA mea- 
surements. 

Components 5 and 4 are modeled analytically as homo- 
geneous synchrotron sources, but when we fit their six- 
frequency spectra using the Stokes / data alone, we found 
that the spectral turnover frequencies, Vpi,ak were not well con- 
strained at low frequency (see the first line in Table 3 for both 
components 5 and 4). However, the fractional linear polariza- 
tion of both components decreased with frequency in a man- 
ner consistent with turnover frequencies closer to 9.0 GHz for 
both components. To quantify this observation we developed 
an analytical model of how fractional linear polarization, m/, 
should decrease with increasing optical depth, including the 
possible effects of depolarization from internal Faraday rota- 
tion. This model is described in detail in the Appendix. 

We repeated the analytical fits including the fractional lin- 
ear polarization data, both with and without the possibility 
of depolarization being due to the observed Faraday rotation 
being internal to the jet. The results for these more highly 
constrained fits are included in Table 3 and plotted in Fig- 
ure |7] panels (a) and (b). To quantify the degree of ob- 
served Faraday rotation in these components we fit a regres- 
sion of the standard form Xohs = RMX^ + Xemit where we found 
RM = -493 ± 23 md m'^ and RM = -250 ± 25 md nT^ for 
components 5 and 4 respectively. These regressions are plot- 
ted in Figure|7] panel (c). Because the total rotations in these 
two cases are relatively small, a regression works well re- 
gardless of whether the observed rotation is internal or ex- 
ternal to the jet. For both components 5 and 4, our results, 
including or excluding the possibility of internal Faraday ro- 
tation, are consistent with peak frequencies of Vpi,ak = 9.0± 1 .0 
GHz and spectral indices of a = -0.5 ±0.1 (we use the con- 
vention S cx i/^"). Because our analysis is not particularly sen- 
sitive to the precise location of the turnover, we adopt these 
values in our subsequent analyses (they are also used to plot 
the spectra which appear in Figure |6]l. 

Component D has an inverted spectrum and is not well 
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Fig. 7. — Stokes / (panel (a)), fractional hnear polarization (panel (b)), 
and electric vector position angles (panel (c)) of components 5 and 4 plotted 
against two homogeneous spectrum models discussed in section 3.3. The 
dashed lines in panels (a) and (b) represent a joint fit to the Stokes / and 
fractional linear polarization data assuming none of the observed Faraday 
rotation is internal to the jet. The solid line in these panels includes a modified 
Bum-style depolarization in the fit assuming that all of the observed Faraday 
rotation is internal to the jet. A A- regression to the observed Faraday rotation 
is plotted in panel (c) as a solid line. 
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TABLE 3 

Spectral Turnover Fits to Components 5 and 4. 



Component 


Method 


^peak 


^peak 


a 


ID 




(GHz) 


(Jy) 




(1) 


(2) 


(3) 


(4) 


(5) 



Stokes I Only 
Stokes I + mi 

Stokes I + mi + FR 
Stokes I Only 
Stokes I + mi 

Stokes I + mi + FR 



Q 511+0.63 
Q A1+0.77 

10.21^.{I 

^■^^-0.69 
^■^^-0.69 



3.94^, 
3.93+^ 

3.06^^' 
3.00+J;' 
3.0l!o: 



"■■'"-0.20 

"■■^■^-0.09 
_n4f;+0.10 

-0.62^1 
-0.46t|l 



Note. — Columns are as follows: (1) Component identifier; (2) Method of fitting 
spectral peak, as described in §3.3; (3) Observed peak frequency in GHz; (4) Peak 
Stokes-I flux density in Jy; (5) Spectral index. 

fit by a homogeneous source spectrum. It can be most 
simply fit by a straight power-law with spectral index a = 
+0.91 ±0.08. In ^3.4.21 we will model component D as an 
inhomogeneous conical jet with power-law dependence on 
magnetic field strength and par ticle density a s described by 
iBlandford & Konigf([T979h and iKonigll (119811) . 

3.3.1. Estimating Magnetic Field Strength 

We can estimate the magnetic field strengths in compo- 
nents 5 and 4 by treating them as homogeneous spheres 
(e.g. Marscher 1987). For a homogeneous volume of plasma 
(iHoman & Wardle..2000. eqn. Al) 



B = 



(l+z) 



(1) 



Where C3 = 5.30 x 10"^ for a = -0.5, fl in mas^, 1^ in GHz, 
S in Jy, and B in Gauss. For an homogeneous sphere at the 
turnover frequency with a = -0.5, r„, = 0.48, Sm is 1.19 times 
the observed flux So, and = (tt/6)9^ where 6d is the angular 
diameter of the component (iHoman & Wardlell2000) . We use 
the approximation, 9d ~ 1.8(00^;,)"'/^, to convert Gaussian 
fitted FWHM dimensions to spherical diameters (iMarscheii 
[T987h . 

Taking the median angular sizes fitted amongst our models, 
we estimate B ^ 0.6 mG for component 5 and B ^ 4.4 mG for 
component 4. These estimates depend on the assumed geom- 
etry and are probably only good to a factor of a few; however, 
our results are not strongly dependent on these values. 

For the inhomogeneous component D, we can apply the 
above technique to make an order of magnitude estimate for 
its magnetic field strength; however, in that case, it would be 
more appropriate to take an optical depth r between 1 and 2, 
and we estimate r ^ 1.5. Using the observed values at 22 
GHz, we estimate B ^ 40 mG. In the following section we 
will allow for a range of possible magnetic field strengths for 
component D. 

3.4. Modeling the Polarization Spectra 
3.4. 1 . Components 5 and 4 

The Stokes / spectra of components 5 and 4 are adequately 
modeled by homogeneous components. Both features show 
strong linear polarization and modest amounts of Faraday ro- 
tation: ~ -35° to -40° and ~ -15° to -20° of total rotation 
by 8.0 GHz for components 5 and 4 respectively. We have no 
direct evidence for how much if any of this observed rotation 
is internal to the jet as neither of these rotations are enough 
to generate large amounts of internal depolarization. In the 



following analysis, we will construct two plausible models, 
one where all of the observed rotation is internal and the other 
where none of the rotation is internal to the jet. 

Our observations indicate that very little, if any, of the ob- 
served circular polarization at 22 and 24 GHz is produced in 
these two features; however, by 8.0 and 8.8 GHz we cannot 
separate the circular polarization between D, 5, and 4. Indeed, 
if we assign all the circular polarization observed at 8.0 and 
8.8 GHz to component D alone, we get very large levels of 
1 .99 ± 0.38% and 1 .55 ± 0.26% respectively. These amounts 
are far larger than the amounts observed on D at 15, 22, and 
24 GHz, where we find less than 1%. So, either component 
D must produce a sharp rise in circular polarization at these 
two frequencies, or components 5 and 4 make a significant 
contribution to the total circular polarization at 8.0 and 8.8 
GHz but not at 15 GHz and above. To accomplish this sec- 
ond scenario, components 5 and 4 must each produce of order 
0.5- 1.0% circular polarization at 8.0 GHz while producing 
< 0.2% at 22 and 24 GHz. 

In reference to the major models of circular polarization 
production considered in 33.11 only models (3) and (5) are 
plausible for these components. Model (2), intrinsic circu- 
lar polarization, needs stronger magnetic field strengths and 
a higher degree of magnetic field order Additionally, the 
spectral dependence of intrinsic CP, ntc oc i^""^, is not steep 
enough. Model (4) requires very large rotation depths which 
are not observed here. 

Model (3) from 33. H is Faraday conversion driven by inter- 
nal Faraday rotation, so for this model, we will assume that 
all of the observed rotation is internal to the jet. In this model, 
the magnetic field order which creates the linear polarization 
at the back of the jet is the same field order which does the 
conversion to circular polarization at the front of the jet. Fara- 
day rotation within the jet rotates the linear polarization by 
angle to allow conversion by this same field, where the cir- 
cular polarization produced is proportional to sin2(/). So for 
this model to operate effectively, we need two components of 
organized magnetic field. The first component is in the plane 
of the sky and first produces the linear polarization and later 
converts it into circular polarization. The second component 
is vector-ordered magnetic field which produces the internal 
Faraday rotation. 

For this model we will assume the first component of mag- 
netic field is provided by a transverse shock. In component 
5, the unrotated electric vector position angle (EVPA) is ap- 
proximately 54°, which gives a magnetic field position angle 
of -36°, nearly perpendicular to the structural position angle 
of -121° for component 5, indeed suggesting that the domi- 
nant field order is due to a transverse shock. For component 
4, the unrotated EVPA is approximately 2°, giving a mag- 
netic field position angle of -88° which is neither perpen- 
dicular nor parallel to its structural position angle of -143°, 
suggesting that the field order in component 4 is due to an 
oblique shock. Our radiative transfer simulation is set-up to 
allow transverse shocks, see 33.2. U however, for the purposes 
of simulating the effects of Faraday conversion in this model, 
the differences between an oblique and transverse shock are 
unimportant, as we only need the shock to provide field order 
in the plane of the sky to first produce the linear polarization 
and then to convert it into circular polarization. 

For the second component of magnetic field which gener- 
ates the internal Faraday rotation, we allow a vector-ordered 
field along the jet axis parameterized as /„ = 0.05 as described 
in 33.2.11 (note that /, = here so there is no contribution 
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from a toroidal field component). With this degree of vector- 
ordered field along the jet axis, the shocks in components 5 
and 4 must be of strength k = 0.64 and k = 0.27 respectively to 
give the observed amounts of linear polarization at high fre- 
quency. W e then put in the magnetic field strengths estimated 
in section ^3.3.11 and adjust the lower cutoff in the relativistic 
particle spectrum, 7/, to give the observed amount of Fara- 
day rotation: 7/ = 19 and 7/ = 8.5 for components 5 and 4 
respectively. The optical depth at each frequency is set by 
the Stokes / spectra determined in ^3.31 For each component, 
we assume a lepton number of £ = 1 .0, representing a pure 
electron-proton plasma. With these parameters set, we com- 
pute the radiative transfer along a line of sight with = 10^ 
cells to minimize stochastic effects from the random compo- 
nents of the magnetic field. The results of this calculation 
appear in Figure [8] 

This model produces the correct sign, amount, and spec- 
trum of circular polarization that we need for components 5 
and 4, with rric in the range of +0.5 to +1.0% at 8.0 GHz and 
falling off to negligible values at high frequency. An impor- 
tant question is: how sensitive are these values to the parame- 
ters chosen in the paragraph above? In short, these values are 
sensitive only to one parameter: the amount of internal Fara- 
day rotation. In the above model we have assumed that all of 
the observed Faraday rotation is internal to the jet. For mod- 
est optical depths, the circular polarization produced by Fara- 
day conversion driven by Faraday rotation is proportional to 
the degree of magnetic field order squared, which is fixed by 
the observations of the linear polarization, and to the product 
TfTc, where tf = Cv''' ™d ''c = CqT the Faraday rotation 
and conversion depth s respectively (Wardle & Homan 20031 
iJones & 0'Dell|[T977l) . The optical depth, r, is also fixed by 
observation, so it is only the coefficients which matter. For 
a = 0.5, the coefficients are given by 



0.2 



1^B± 7/ 



cot 61' 



and 



7/ 



(2) 



(3) 



where 7 = y/ v/vb± is the 7 of the radiating particles, and 
and are constants of order unity. Note that the de- 
pendence on /„ given for is approximate and only good for 
small values of /„ < 0.1. 

In summary, if all of the observed Faraday rotation is as- 
sumed to be internal, then rp is also fixed by observation 
and the combination of parameters £f„B~^ ln7;/7,^ is likewise 
fixed. Thus, the only place where parameter choice affects the 
results of this model is through the conversion depth tc which 
depends only weakly on the magnetic field strength, B± and 
the lower cutoff, 7/, through the logarithm In (7/7;). In this 
respect, the rotation driven conversion model is attractive as 
the sign, amount, and spectrum of circular polarization pro- 
duced is nearly fixed by the observed properties of the linear 
polarization. 

In contrast, the helical magnetic field model, model (5) 
from ^3.11 has fewer constraints because here we assume that 
all of the observed Faraday rotation is external to the jet, so as 
to assess the ability of the helical field itself to produce the ob- 
served circular polarization. As above, we assume /„ = 0.05 
for the vector-ordered component of magnetic field along the 
jet axis, and with this value, we require a toroidal field of 
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Fig. 8. — Linear and circular polarization spectra for comp onents 5 and 4 
from the rotation driven conversion model described in j|3.4| Panel (a) plots 
the fractional linear polarization model against the data, panel (b) plots the 
electric vector position angle (EVPA) against the data, and panel (c) plots 
the predicted fractional circular polarization as a function of frequency. The 
model EVPA values in panel (b) have been shifted by a constant offset to 
align with the high frequency data. Note that in panel (b), the EVPA for 
component 4 has been shifted by —120 degrees to allow it to fit comfortably 
on this plot. 
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/, = 0.56 and /, = 0.78 to match the observed levels of linear 
polarization in components 5 and 4 respectively. We note that 
the unrotated EVPA of component 4, computed above, is not 
consistent with a toroidal field which should be perpendicular 
to the jet axis; however, for the purpose of this analysis we 
will ignore this inconsistency. 

With /, 3> /«, the overall field-order is a high pitch-angle 
helix as sues ested bv lGabuzdaetaLl (|2008|) for 3C279, and 
we can obtain similar levels of circular polarization as the 
rotation driven conversion model by choosing 7; = 20 and 
7/ = 10 for components 5 and 4 respectively. We have also 
taken I = 10"^ to completely eliminate any internal Faraday 
rotation, and no shock has been assumed for either compo- 
nent, k= 1.0. For a helical field, where the field order varies 
across the jet cross-section, a single line of sight is not suf- 
ficient to determine the emergent polarization, so we take a 
100 X 100 X 100 cube and set to zero any cells outside a cylin- 
drical jet cross-section. The emerging polarization for both 
components is plotted in Figure |9] 

As with the rotation driven conversion model, the helical 
field is able to produce circular polarization with the correct 
amplitude, sign, and spectra; however, the amount and sign of 
the produced circular polarization are no longer determined 
by the linear polarization and are the result of choices we 
made about field direction and the magnitude of 7;. The spec- 
trum of the circular polarization from the helical field falls off 
more slowly than for the rotation driven conversion case, but 
the levels are small enough at 15 GHz and above to still be 
acceptable. This shallower spectrum is due to the fact that 
the helical field produces CP via pure Faraday conversion and 
depends only upon tc, whereas rotation driven conversion de- 
pends on the product of t/tTc which includes two powers of 
the optical depth. 
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3.4.2. Component D 

Component D represents the base of the jet as observed 
in our VLBA images, and it has an inhomogeneous spec- 
trum which is well fit as a single power-law S oc v^" where 
a = +0.91 ± 0.08 as described in 33.31 Its linear polarization 
is approximately flat or perhaps slightly decreasing with fre- 
quency from nil — 3.6 to 2.3% from low to high frequency. 
The EVPA of the linear polarization stays roughly constant 
between +110° to +120° except for 24 GHz, where we find 
+101° ± 3°. Its circular polarization at 22 and 24 GHz is 
0.86 ± 0.20% and 0.94 ± 0.22% respectively, and if we assign 
all the observed circular polarization at 15 GHz to component 
D, we find a similar amount of 0.88 ± 0.15%. 

These properties are all suggestive o f an inhomogeneous 
component of the type i nvestigated by iBlandford & Koni"gl 
([19790 and lKonigll(ll981h where the magnetic field and parti- 
cle densities scale as power laws in a conical jet with Z?x r~"' 
and K cx r~". The EVPA of the linear polarization described 
above corresponds to a jet magnetic field which is neither per- 
pendicular nor parallel to the jet axis. The later is assumed 
to lie somewhere between the structural position angles of 
components 5 and 4 of -121° and -143° respectively. Ad- 
ditionally, the observed levels of linear polarization at high 
frequency (2-3%) are too low to generate the nearly 1% cor- 
responding circular polarization if the linear polarization is 
taken as a direct measure of the field order. In all likelihood, 
Faraday rotation and depolarization, either internal or exter- 
nal to the jet, are responsible for the offset EVPAs and the 
low observed levels of linear polarization. If the culprit Fara- 
day screen is external to the jet, then it must scale in a similar 
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Fig. 9. — Linear and cii'cular polaiization spectra fo r com ponents 5 and 
4 from the helical field conversion model described in j|3.4| Panel (a) plots 
the fractional lineal' polarization model against the data, panel (b) plots the 
electric vector position angle (EVPA) against the data, and panel (c) plots 
the predicted fractional circular polarization as a function of frequency. The 
model EVPA values in panel (b) have been shifted by the same constant offset 
used in Figure [8] Note that in panel (b), the helical field model predicts a 
flat EVPA with frequency because none of the observed Faraday rotation is 
assumed to be internal to the jet. As in Figure [8] the EVPA for component 4 
has been shifted by —120 degrees to allow it to fit comfortably on this plot. 
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power-law fashion in order for the jet properties themselves 
to maintain an approximately constant level of depolarization 
and rotation with frequency. While such a scaling of an ex- 
ternal screen is certainly possible, it seems more likely that 
the depolarization and rotation is occurring internal to the jet, 
where it naturally would scale in the appropriate manner with 
B and K. 

In constructing our inhomogeneous model of component D, 
we will assume that all of the observed depolarization and ro- 
tation are occurring internal to the jet, matching our assump- 
tion for the rotation-driven conversion model of components 
5 and 4. This assumption provides the most stringent con- 
straints on component D, as the same magnetic field model 
must produce the correct signs and amounts of circular po- 
larization, Faraday rotation, and depolarization. From 33.11 
models (2), (3) and (4) for circular polarization production 
all fit this scenario, and we can simultaneously investigate 
all three by constructing inhomogeneous jets consisting of 
vector-ordered field along the axis and disordered field which 
may be shocked from unit length to length k as described in 

We simulate conical, inhomogeneous jet emission by run- 
ning our radiative transfer simulation for many lines of sight 
through the center of a jet where each line of sight corre- 
sponds to a new radius, r, from the base of the jet. The mag- 
netic field strength and frequency where r = 1 are scaled at 
each radius, r, in a power-law fashion given by Konigl ( 1981) 
where B± cx r"'" and v cx r"*'", where k„ is given bv Konigli 
(Il981h and corresponds to a power-law scaling of not only 
B± but also the particle density, K (x r~". In all of our models, 
we assume an optically thin spectral index of a = -0.5. For 
m = 1 .8 and « = 2.8, =1.8 and we find a good matc h to the 
Stokes / spectrum from our model. For these values, Konigli 
(Il981h would predict S cx iy"*""^^, which agrees well with our 
simulation. Other combinations of m and n could also match 
our Stokes / spectrum, and we explore three such combina- 
tions. 

We set the magnetic field strength and degree of vector- 
ordered field at an observed frequency of 22 GHz. We require 
the vector ordered field to scale like to conserve magnetic 
flux, and we adjust the scaling of the random component of 
the field to give the overall correct scaling for B± cx r~"\ 
For each location in the jet we compute the emerging radi- 
ation at our VLBA observing frequencies and integrate the 
results from the entire inhomogeneous component, cutoff at 
radii which produce low and high frequencies of 1 GHz and 
100 GHz respectively in the observer frame. The Stokes / 
spectrum is then scaled by a single multiplicative factor, com- 
mon to all frequencies, to best align with our observations. 
The model EVPA is also rotated by the same angle at all fre- 
quencies to best align with the observed EVPAs between 15 
and 24 GHz; however, this angle is not arbitrary as it tells us 
the direction of the net magnetic field in the jet of 3C279 and 
should correspond to a sensible value for the magnetic field 
model of the simulation. The final result is a model spectra 
of Stokes /, fractional linear polarization, EVPA, and circular 
polarization. 

We initially ran a coarse grid of models with m = 1.8 and 
n = 2.8, exhausting all combinations of the following parame- 
ter values: £= 1.0,0.1,0.01, /„ = 0.1 to 0.9 in 0.1 steps, B± = 
0.1,0.05,0.02,0.01 G, 7, = 3,5,8,10,15,25,50, and;t = 0.1 
to 1.0 in 0.1 steps. We judged degree of agreement between 
model and data by comparing the fractional circular polariza- 
tion, fractional linear polarization, and EVPA at 15, 22, and 



24 GHz. We don't expect every variation and wiggle in the 
observed data to be reproduced by the inhomogeneous power- 
law models. In fact, we expect the inhomogeneous models 
to produce smoother-than-observed values for the polariza- 
tion; however, we should be able to reproduce the general 
trends and levels observed. In this light, we required "plau- 
sible" models to produce an average circular polarization in 
the range = 0.8 to 1.0%, with no value smaller than 0.6% 
or larger than 1 .3%, and we required average linear polariza- 
tion in the range mi = 2.2 to 2.5%, with no value smaller than 
2.0% or larger than 2.8%. These restrictions will only allow 
models to be plausible if they reproduce the general levels and 
slopes in the observed data. Out of more than 7000 models in 
this coarse grid, we found only two that fit these criteria; both 
models had £ = 1.0, B± = 0.1 G, /„ = 0.7-0.8, 7/ = 5, and 
/t = 0.8-0.9. 

Physically, these plausible models represent a jet dominated 
by vector-ordered magnetic flux along the jet axis, produc- 
ing large amounts of internal Faraday rotation/depolarization 
with a significant contribution from intrinsic circular polariza- 
tion as well as rotation-driven conversion. A range of similar 
models, while not meeting our plausible criteria did show sim- 
ilar trends in mi and mc, and we wished to explore this type of 
physical model with a finer grid. This new grid looked at ev- 
ery combination of ^ = 1 .0, 0.5 , 0. 1 , B± = 0. 1 , 0.05 , 0.02, 0.01 
G, /„ = 0.50 to 0.95 in 0.05 steps, and 7/ = 3,4,5,6,7,8,9, 10. 
We realized that the shock strength, k, was not a physically 
meaningful parameter in this case, so we set A: = 1.0 for our 
fine grid. Out of these 900+ models, we found 4 models meet- 
ing the plausible criteria discussed above. All of these models 
had e=l.O,Bj_ = 0.05G or O.IG, /„ = 0.7-0.9, and 7, = 5 or 
6. The best overall model from this set is plotted in figure [TOl 

The model plotted in Figure [TO] has ^ = 1.0, B_l = 0.05 G, 
/j, = 0.8, and 7; = 6. While the comparison between model 
and data was done at 15 GHz and above, we show here the 
full spectra of the model at all our observing frequencies. 
Note that the model deviates from the fractional linear po- 
larization and EVPA data below 15 GHz. This trend could 
be easily explained if the amount of internal rotation reduced 
more quickly than expected at larger r, perhaps due to local 
conditions in the jet or an increase in the lower cutoff to the 
relativistic particle spectrum, 7/, with r. This deviation illus- 
trates an important limitation of this model which can only 
simulate broad spectral trends from smoothly varying physi- 
cal properties. 

For the circular polarization at 15 GHz and below in Figure 
[TOl we have added a dotted line which includes the contri- 
butions from the rotation-driven conversion models of com- 
ponents 5 and 4, given in figure [8] Note that the total cir- 
cular polarization at low frequency is approximately correct 
with this addition, although the total falls a bit higher than the 
measured 8.0 GHz value. 

The intrinsic magnetic field direction in this model corre- 
sponds to a position angle of -129°, which places it between 
the structural position angles of components 5 and 4. This 
is consistent with a parallel magnetic field ordered along a 
jet axis which points at ~ -129°. All of our plausible mod- 
els gave de-rotated magnetic field directions between -129° 
and -133° to give the best match to the 15, 22, and 24 GHz 
data. Given that all of these models are dominated by vector- 
ordered field along the jet axis, i.e. /„ > 0.7, this is an im- 
portant consistency check, revealing that not only does this 
model produce the appropriate amounts of both linear and cir- 
cular polarization, but it also produces approximately the right 
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Fig. 10. — Full po larization spectra of the inhomogeneou.s component D. The data are plotted a.s open squares against a dashed line representing our best model 
described in j|3.4.2l As described in that section the quality of the agreement was judged only at 15, 22, and 24 GHz; however, we plot the full spectrum here. 
Panel (a) is Stokes /. Panel (b) is fractional circular polarization. In panel (b), the 8.0 and 8.8 GHz values ai'e plotted as soHd squai'es to indicate that they may 
include significant contributions from components 5 and 4. The results from our rotation driven conversion model in figure[8]are added to the inhomogeneous 
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amount and sign of EVPA rotation for the emerging radiation. 

To check our dependence on the power-law exponents, m 
and n, we repeated the coarse grid described above for two 
additional combinations, m = 1.3, n = 3.4 and m = 2.3, n = 
2.3, both of which match our Stokes / spectrum. We found 
only one model which fit our plausibility criteria, and it had 
m = 2.3, n = 2.3, £ = 1.0, B_l = 0.02 G, /„ = 0.8, 7; = 10, 
and k = 0.7. The intrinsic magnetic field direction in this 
model corresponds to a position angle of -128°. These re- 
sults present a very similar physical picture to the best models 
with OT = 1 .8 and n = 2.8. 

4. DISCUSSION 

4.1. Homogeneous Components 5 and 4 

For components 5 and 4, we find little to no circular po- 
larization at 22 and 24 GHz and as mu ch as 0.5% to 1.0% 
circular polarization by 8.0 GHz. In ^3.4. II we explored 
two possible physical models to explain this circular polar- 
ization spectra both based on Faraday conversion of linear 
polarization into circular The first model corresponds to 
model (3) described in section 33.11 rotation-driven conver- 
sion, and involves linear polarization generated by a shocked 
magnetic field, a vector-ordered component of field along the 
jet axis which causes Faraday rotation of the linear polariza- 
tion within the jet, and finally conversion of the rotated linear 
polarization into circular polarization by the shocked mag- 
netic field. The results of this model were plotted in Figure 
[8] The second model corresponds to model (5) described in 
section 33.11 conversion in a helical field, and involves linear 
polarization generated at the back of the jet being converted 
into circular polarization by the magnetic field at the front of 
the jet which is at some angle with respect to the field at the 
back of the jet. 

Both of the models for components 5 and 4 described above 
are Faraday conversion models, as this will give the neces- 
sary steep spectrum, and both explain the circular polariza- 
tion of these components adequately. However, we prefer the 
rotation-driven conversion model because it not only explains 
the circular polarization but also the linear polarization, in- 
cluding the amou nt and sign of the observed Faraday rotation. 
As described in 33.4.11 in this model the amount and sign of 
the observed circular polarization is essentially fixed by the 
amount of internal Faraday rotation along with the observed 
degree of magnetic field order given by the linear polarization. 
On the other hand, the helical field model has greater freedom 
to tune the amount and sign of the observed circular polariza- 
tion because it decouples the observed Faraday rotation into 
an external screen which may have different properties than 
the jet itself. 

If all of the observed Faraday rotation in components 5 
and 4 is internal to the jet, the Faraday conversion depth at 
8.0 GHz, Tc, is constrained to lie in the ranges 1.5-3.0 and 
0.75- 1.5 for components 5 and 4 respectively to generate 
fractional circular polarization in the range mc = 0.5- 1.0% 
for each component. For the observed optical depths of 
these components, the Faraday conversion depth constrains 
the Lorentz factor ratio between the emitting particles and the 
low energy cutoff, 7/7;, to be in the range 20-400 for com- 
ponent 5 and in the range 4.5-20 for component 4. Using 
the estimated magnetic field strengths from 33.3. ll to find 7 
at 8.0 GHz, the lower cutoff in the particle energy spectrum 
is constrained to be 1 < 7/ < 27 and 9 < 7/ < 43 for compo- 
nents 5 and 4 respectively. Note that these limits scale with 



the estimated magnetic field strength like 1 / \/B. 

From the observed Faraday rotation itself, we obtain only a 
joint constraint on futlnji/'^^ as described in 33.4.11 If there 
is significant cold "thermal" matter inside the jet, this con- 
straint will be modified by an additional term: +/?f(3/4)/7/. 
where R,- is the number ratio of cold to relativistic particles in 
the jet (iJones & 0'DellllT977l) . 

4.2. Inhomogeneous Component D 

The inhomogeneous component D lies at the very base of 
the parsec scale radio jet, and we argue in 33.4.21 that its 
S tokes / and polariza t ion sp ectra can be well m odeled by 
a iBlandford & Konigll (Il979h and iKonigll (1198 ll) style con- 
ical jet with magnetic field and particle strengths that scale 
as power-laws with distance from the jet origin. We prefer a 
physical model for component D where the linear polarization 
is depolarized and rotated by the same magnetic field and par- 
ticles in the jet which produce the Stokes / emission, linear 
polarization, and circular polarization. These circumstances 
provide the strongest constraints on the physical properties of 
the jet magnetic fields and particles as these properties alone 
must produce consistent amounts and directions of both lin- 
ear and circular polarization across multiple frequencies. The 
physical model required by these constraints is dominated by 
a vector-ordered magnetic field along the jet axis with /„ > 0.7 
which does three things: (1) produces the observed linear po- 
larization, (2) generates internal Faraday rotation which ro- 
tates the EVPA and depolarizes the linear polarization to the 
low fractional levels observed, and (3) produces the observed 
circular polarization largely as intrinsic circular polarization, 
although there is some contribution from rotation-driven con- 
version of the linear polarization into circular. 

It is interesting to note that the dominant mechanism for 
producing circular polarization in components 5 and 4 is ro- 
tation driven conversion, not intrinsic circular polarization; 
however, these homogeneous components are not dominated 
by vector-ordered magnetic field along the jet axis which 
would be required for efficient intrinsic circular polarization 
production. In fact, for components 5 and 4, some fraction 
of vector-ordered field, /„, is required to generate the neces- 
sary internal Faraday rotation, and indeed the same polarity 
of vector-ordered field is required in components 5 and 4 as 
we deduce for component D. 

4.2.1. Estimating the Jet Magnetic Flux 

For our model of the inhomogeneous component D, we re- 
quired the vector-ordered field to scale as r~^ along the jet 
axis, consistent with conservation of magnetic flux in a coni- 
cal jet. Continued conservation of magnetic flux further down 
the jet in components 5 and 4 would leave some smaller frac- 
tion of vector-ordered field in these components with the same 
polarity. For components 5 and 4, it is not possible to com- 
pute the required magnetic flux because, as described in 34.11 
/„ is tied to other unknown parameters, and it is also not clear 
whether components 5 and 4, which are on different structural 
position angles, occupy the entire jet width at their location or 
represent bright spots within a jet with a wider opening angle. 
However, these limitations don't apply to component D which 
has /„ > 0.7, and we can estimate the jet magnetic flux by fit- 
ting the transverse size of component D in our VLBA data. 
We find the transverse size to be in the range 0.02-0.04 mas; 
however, we note that we are at the very limits of our resolu- 
tion here and can only use this value to estimate the jet mag- 
netic flux. To convert from Gaussian FWHM to diameter of 
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the jet, we multiply by 1 .8 (lMarscherlll987h . and at a redshift 
z = 0.536, 1 mas corresponds to a linear scale of 6.3 parsecs 
for standard cosmological parameters of Hq = 70 km/s/Mpc, 
17m = 0.3 and f^A = 0.7. For /„ = 0.7 and B±_ = 0.050 G, the net 
magnetic flux in the jet is in the range ~ 2 x 10-''*- 1 x lO^*^ G 

2 

cm . 

This net magnetic flux is, to our knowledge, the first esti- 
mate of the magnetic flux carried by a jet that takes proper 
account of the vector order of the magnetic field. In the ab- 
sence of field line reconnection or entrainment, it is a con- 
served quantity and is therefore equal to the net poloidal mag- 
netic flux at the central engine. If the field originates in 
the black hole magnetosphere, then the expected magnetic 
flux for an Eddington black hole is 2 x cm'' 
(Wardle & Homan 2003, after correcting an arithmetic er- 
ror; Begelman, Blandford & Rees 1984). Setting this equal 
to the measured magnetic flux implies a black hole mass in 
the range - 2 x 10''-6 x lO'' M© in 3C279. This would be 
at the top of the range of me asured or inf erred black hole 
masses found in the literature dWoo & Urrvl l2002) but is not 
implausible; however, this estimate is an order of magni- 
tude larger than the inferred black hole mass for 3C 27 9 itself 
(IWoo & Urrvll2002HWang et"ani200linuetani2006h . Alter- 
natively, the poloidal magnetic field may be anchored in the 
accretion disk, generally leading to larger magnetic fluxes be- 
cause the area is much la rger (e.g. iBlandford & Pavnel[T982t 
IVlahakis & Koniglll2004 . In the latter paper, the net mag- 
netic fluxes in the two models presented are lO"'^ G cm^ and 
lO-'^ G cm- (Vlahakis, private communication) which span 
our estimated range. The concordance between estimated 
and expected values suggests that circular polarization mea- 
surements can indeed measure a fundamental property of the 
central engine, and may even enable discrimination between 
models. 

4.2.2. Jet Composition 

The other physical parameter of interest constrained by this 
model is the lepton number I, which we find to be 1 .0 in all 
of our plausible models; however, the next smallest value we 
tested was I = 0.5. The lepton number is most strongly con- 
strained by the circular polarization which is largely due to 
instrinsic circular polarization in this model. Smaller val- 
ues of £ require larger values of B, such that B cx l'^, so if 
£ = 1 .0 for B ~ 0.05 G as we find here, I = 0.5 would corre- 
spond to B ~ 0.2 G.^ Pushing this further, values as small as 
£ = 0.1 might be allowed for stronger magnetic fields up to a 
few Gauss; however, larger values of B produce correspond- 
ingly larger values of the magnetic flux which is already quite 
large when estimated assuming B = 0.05 G, see 34.2.11 Ad- 
ditionally, from our observed flux density, angular size, and 
estimated optical depth at 22 GHz, we estimated B ~ 0.04 G 
in 33.3.11 and while differing from this value by a factor of 
a few is certainly possible, the large deviations required for 
l = Q.\ or less are implausible. We therefore find that in this 
model, the jet in 3C279 is at least dynamically dominated by 
protons with I likely to be greater than or equal to about 0.5, 
corresponding to > 75% of the radiating particles being elec- 
trons. 

In the first parsec-scale circular polarization study of 

* Note that the dependence of B on is further compHcated by the fraction 
of uniform field, /„, such that the quantity fulB'^-^ is approximately constant 
for /j, near unity; however, as /„ is already > 0.7 not much additional reduc- 
tion in I can be gained by increasing f\,. 



3C 279, IWardle et all (Il998h preferred a jet dominated by 
electron-positron pairs, and their argument was based on their 
conclusion that 7; < 20 and therefore a jet dominated by 
electron-positron pairs was required by ki netic luminosity ar- 
guments originated by Celotti & Fabian (1993). Applying 
the same calculation as Wardle et al. (1998) for jet energy 
flux to these observations, we derive a similar constraint with 
Fe = 1.33 X 10^^[l + 153£/7/] ergs/s. This calculation assumes 
a magnetic field st rength of B± = 0.05 G and, like that of 
IWardle et al.l ( 1199 8'). assumes equipartition between the mag- 
netic field and particle energies. Therefore, this calculation 
should be taken as a lower limit on the jet kinetic luminosity 
when the field sti-ength is 0.05 G. IWardle et al.l (11998.) found a 
value for Fe about twice this, and they concluded that 7; < 20 
required a mostly electron-positronjet to avoid carrying much 
more energy than appears to be dissipated on larger scales. 
It should be noted that in our model here, because the field 
strength, B, scales roughly with i'^, reducing I far below 1 .0 
does not reduce the energy carried by the jet, but rather in- 
creases it sharply as t~^. 

It seems that in this picture for component D there are only 
two ways to reduce the energy carried by the jet in the cal- 
culation given above. The first is reducing the magnetic field 
strength, which is possible as we had at least one model with 
B±^ = 0.02 G when we tried power-law indices of m = 2.3 and 
n = 2.3, as described in §3.4.2. The second i s by increasing 7/. 
As we not ed in the introduction, wo r k bv [ Beckert & Falckd 
(I2OO2I) and lRuszkowski & BegelmanI (l2002b showed that the 
addition of thermal plasma to the jet generates enough internal 
Faraday rotation to allow larger values of 7/ therefore relax- 
ing the constraint on 7/ found by Wardle et al. ( 1998). In the 
model for component D explored here, a similar ambiguity 
exists for the value of 7/. All of our plausible models have 
7/ < 10; however, this constraint on 7; is driven by the need to 
generate enough internal Faraday rotation and depolarization 
by the relativistic plasma itself and the addition of thermal 
matter to the jet would relax this constraint on 7/. Our models 
for components 5 and 4 provide separate limits of 5 < 7; < 35, 
based on Faraday conversion, that are not subject to this ambi- 
guity. Taken together, the uncertainties in the value of B±^ and 
7/ allow for smaller values of Fe, p erhaps an order of m agni- 
tude smaller than that calculated bv IWardle et al.l (Il998l) . 

4.2.3. Uniqueness of Our Model 

The model described for component D in the previous sec- 
tions is of an inhomogeneous jet which is dominated by 
vector-ordered, poloidal field along the jet axis. In this model, 
the poloidal magnetic field is responsible for generating the 
observed circular polarization as intrinsic circular polariza- 
tion, and the same field provides the internal Faraday rotation 
and depolarization which produces the relatively low levels 
of observed linear polarization with a distinctly rotated EVPA 
which is offset from the jet axis. The inhomogeneous nature 
of the component explains the Stokes / spectrum, the circu- 
lar polarization spectrum, and the approximately flat spectrum 
for the fractional linear polarization and its EVPA, despite the 
large Faraday depth required. It is these properties that we 
consider to be robust aspects of our model, and they lead nat- 
urally to our constraints on the composition of the radiating 
particles as parameterized by lepton number £. 

We expect that the details of the inhomogeneous model we 
have used are not unique, and that other inhomogeneous mod- 
els could reproduce our observations just as well if they have 
the essential properties listed above. Indeed, as discussed 
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in §3.4.2, any inhomogeneous model with smoothly varying 
magnetic field and particle densities is itself an idealization of 
the base of the jet which likely has local variations in these 
properties. In this respect, we see our detailed numerical re- 
sults from §3.4.2 only as estimations of the physical proper- 
ties of this feature. 

Finally, we reiterate that if an external Faraday screen is re- 
sponsible for rotating and depolarizing the linear polarization, 
a wider range of general magnetic field and particle models 
could produce the observed circular polarization. Parsec-scale 
Faraday rota tion measure observ atio ns of 3C279 have been 
pubHshedjw lT^v^ (IT9981 llOOOh and lZavala & Tavkii (I20QTI 
I2OOI 12001 spanning five epochs from early 1997 through 
mid-2001. Their results agree with ours in that they find con- 
sistently negative, but variable, Faraday rotations in the core 
region of magnitudes encompassing those we see in compo- 
nents D, 5, and 4. The variable rotation measures they observe 
fit naturally into our model where the amount of poloidal mag- 
netic field drops off sharply from the core region, and newly 
emerging components can sample a range of internal rotation 
measures as they propagate down the jet. However, an exter- 
nal screen, very close to the jet, could also have these proper- 
ties. Further out in the jet, these authors find a strong jet com- 
ponent to have a small amount of negative Faraday rotation at 
about 3 milli-arcseconds from the core in their 1997 and 1998 
epochs; however, by 2000, this component shows a small 
amount of positive Faraday rotation, less than +50 rad/m^ in 
the observer's frame of reference. This shift to positive rota- 
tion measure could be related to a collision and re-ahgnment 
event this component underwent in mid- 1998 (Hom an et al.l 
12003), or the component could simply be far enough from the 
central engine that other external Faraday screens between us 
and the source are playing the dominant role.' 

This evidence indicates that the negative Faraday rotation 
observed near the core region by ourselves and others is in- 
deed produced local to the jet. It could be due to the fields 
internal to the jet, as we propose here. Alternatively, it could 
be generated externally, but close to the jet, by fields and par- 
ticle densities which scale with frequency and position in just 
the right way to produce both the approximately constant frac- 
tional linear polarization and offset EVPA observed in com- 
ponent D at all frequencies. We regard the fact that linear 
polarization of this kind comes out naturally from the internal 
rotation model as strong evidence in favor of it. 

4.3. Comparison to Integrated UMRAO Monitoring 

An interesting result from the integrated UMRAO monitor- 
ing of 3C279 at 4.8, 8.0, and 14.5 GHz is that 3C279 shows 
occasional sign reversals in it s circular polariza tion at 4.8 GHz 
but not at higher frequencies (^ Aller et ani2006 . and Figure 5). 
An important question is whether these sign reversals can be 
understood in the context of the model of the core region pre- 
ferred for the VLBA data. From Figure |6] it is clear that the 
core region flux at 4.8 GHz is dominated by components 5 
and 4. The Faraday conversion models preferred for compo- 
nents 5 and 4 will naturally produce negative circular polar- 
ization at high optical de pths, at a f requency about one-third 
of the spectral turnover (Jones & O'Dell 1977). With spec- 
tral turnovers near 9.0 GHz, components 5 and 4 should not 

' Note that lZavala & Tavloil 120041) find an integrated rotation measure of 
+3 1 rad/m^ for 3C279 which is nearly an order of magnitude smaller than 
even the smallest negative rotation measure we find here for component 4. 
and thus correction for this integrated value would have essentially no impact 
on our results. 



be producing negative circular polarization at 4.8 GHz in the 
VLBA epoch investigated here, and indeed negative circular 
polarization is not observed near this epoch; however, when 
these components were first originating near the base of the 
jet, they must have been more compact with higher turnover 
frequencies and may well have generated the negative circular 
polarization observed at 4.8 GHz in mid-2005 along with the 
increase in Stokes / at that time. Strong negative circular po- 
larization is also seen at 4.8 GHz at several other epochs dur- 
ing this period: early 2004, mid-2004, mid-2006, late 2006, 
and late 2007. Each of these events also appear to be linked 
with increases in Stokes /, suggesting that new components 
may be ejected from the base of the jet at these times. 

5. CONCLUSIONS 

We have made multi-frequency, parsec scale observations 
of the core region of 3C 279 in Stokes /, linear polarization, 
and circular polarization. The core region of 3C 279 consists 
of three main components: two homogeneous components, 5 
and 4, and the inhomogeneous component D at the base of the 
jet, and we have modeled the full polarization spectra of these 
components with radiative transfer simulations to constrain 
the magnetic field and particle properties of the parsec-scale 
jet in 3C279. 

We find that the polarization properties of the core region, 
including the amount of linear polarization, the amount and 
sign of Faraday rotation, and the amount and sign of cir- 
cular polarization can be explained by a consistent phys- 
ical model. The base of the jet, component D, is mod- 
eled as a co nical inhomogeneous Blandford-Konigl style jet 
(iBlandford & Konigl 1979; Konigl 1981). The magnetic field 
of this feature is dominated by vector-ordered poloidal mag- 
netic field along the jet axis, and we estimated the net mag- 
netic flux of this field. This poloidal field at the base of the jet 
produces intrinsic circular polarization and depolarizes and 
rotates the linear polarization. In the homogeneous compo- 
nents 5 and 4, further down the jet, the magnetic field is 
dominated by local shocks, and a much smaller fraction of 
vector-ordered, poloidal field remains along the jet axis. This 
remaining vector-ordered field provides a consistent internal 
Faraday rotation, which allows Faraday conversion to gener- 
ate the appropriate amount and spectra of circular polarization 
from these components. 

We note that this physical model is not unique if one allows 
the observed Faraday rotation and depolarization to occur in 
screens external to the jet. Such external screens would de- 
couple the linear and circular polarization and allow a wider 
range of physical parameters; however, we find the fact that 
all of the essential polarization characteristics can be simul- 
taneously produced by magnetic field and plasma properties 
internal to the jet itself to be compelling motivation for this 
model. 

With this model, we can additionally constrain the particle 
properties of the jet. We find the jet to be kinetically domi- 
nated by protons with a lepton number I > 0.5 corresponding 
to > 75% of the radiating particles being electrons, and there- 
fore we cannot rule out a significant admixture of positrons. 
Based on the amounts of Faraday conversion deduced for the 
homogeneous components 5 and 4, we find a plausible range 
for the lower cutoff in the relativistic particle energy spectrum 
to be 5 < 7, < 35. 
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APPENDIX 

OPACITY EFFECTS ON LINEAR POLARIZATION 

Here we explore the effect of optical depth, t, on the observed fraction of linear polarization, nii, from a partially ordered mag- 
netic field. Our goal is to obtain an analytic model for mi as a function of t which can be used in fitting spectra of homogeneous 
synchrotron sources. We take the ordered part of the magnetic field to be in the plane of the sky at a position angle of 90° so 
that the emitted linear polarization is entirely Stokes-Q. We also assume no internal Farad ay rotation o r conversion. Under these 
circumstances, the equations of radiative transfer reduce to the following expressions fe.g. lJonelll 19881) . 

For a completely ordered magnetic field, we have the exact differential equations ... 

^+I+CqQ = J (A1) 

where eg = 0.6923 and (q = 0.7500 are the emission and absorption coefficients respectively for a = -0.5. For optically thin 
emission in a completely ordered magnetic field, the fractional linear polarization is equal to the emission coefficient: m/(T = 0) = 

For a partially ordered field, we will substitute appropriately scaled emission and absorption coefficients, e and for eg and 
Ce respectively. 

^ + I+CQ = J (A2) 
dr 

dQ 

-^ + Q + CI=eJ 
dr 

The scaled emission coefficient e must give the observed fractional linear polarization at r = 0, so e = m/(T = 0), and we can 
guess that the appropriate scaling for the absorption coefficient is similar. We parameterize C = Ce ( '"'i^'^""' ) ' where 77 is an 
unknown power, likely to be close to unity. 

These coupled, first order differential equations can be easily solved and we find the resulting fractional linear polarization as 
a function of opacity. 

^ ^ e ^ (1 + e)(l - 0(1 - - ( 1 - £)( 1 + C)( 1 - e-^'-O-) 
/ {\ + e){\-0(l-e-''^+Or) + {i-e){\+C){\-e-(^-0^) 

In the limit r ^ 0, m; = e as we expect, and as t ^ inf, mi = jz^. All of these expressions are exact for a completely ordered 
magnetic field; however, we are interested in a partially ordered field using the scaled versions of the emission and absorption 

coefficients given above: e = m/(r = 0) and C = Cg ■ 

To test these expressions, we ran the radiative transfer simulation described in ^3.2l with a = -0.5 for various fractions, /„, for 
the vector-ordered magnetic field, taken to be in the plane of the sky. We used = 10^ cells along the line of sight to minimize any 
net contribution from the random component of magnetic field in each cell. We also set the lepton number £ = 10"^ to eliminate 
the effects of internal Faraday rotation. The fractional polarization, mi from these runs is plotted against optical depth, r, in figure 
lAll The model given in equation lA3l is plotted for each value of /,, where the best values of e and were found by a least-squares 
fit. The best-fit values for e and C, are given in the legend of figure IaTI and we can see that indeed e = m;(r = 0) as expected, and 

we find that our scaled value for C = Ce ^ '"'(t-=0) ^ excellent approximation if 77 = 0.86. 

An additional factor to consider is the possibility of depolarization due to internal Faraday rotation. 'Burn' ('1966') predicts 
additional depolarization of sin(<I>)/<l' where ^Bum = 2.0Ax; however, this is for a purely optically thin case. Figure A2 shows 
the results of our radiative transfer simulation with /„ = 0.3 with the jet oriented at 6* = 45° to the line of sight to generate internal 
Faraday rotation. The amount of Faraday rotation as a function of optical depth was chosen to be similar to what we observe for 
components 5 and 4, with < 40° of internal rotation up to an optical depth of r ^ 1 . Note that the internal rotation is well fit by 
a A- regression up to Ax — 40°. Panel (b) of this figure shows the additional effect of internal depolarization beyond the pure 
optical depth reduction to m/. To include the effect of internal Faraday depolarization in these plots, we used the linear regression 
fit to Ax = RMX^ to produce the $ values. Note that $b„„, somewhat over-predicts the effect of internal depolarization, and we 
found that a modified factor $ = 1 .7 Ax produced a better match to the simulation data. 

Taken together, equation A3 and above depolarization rule for internal Faraday rotation, allow us to include fractional linear 
polarization, mi, in our analytical fit to the observed Stokes / spectrum of homogeneous jet components to find their peak 
frequency. To test this technique, we ran this procedure on a detailed numerical simulation of component 5 from our radiative 
transfer program. We took the results of this simulation at each frequency and treated them as data with the same uncertainties 
on mi and Stokes / as the real data from component 5. The simulation had /„ = 0.05 with the jet oriented at 6* = 1 .5° to the line of 
sight with a bulk Lorentz factor of F = 15. To generate the right amount of fractional linear polarization at high frequency a shock 




Fig . A 1 . — Fractional linear polarization plotted against optical depth from our radiative transfer simulat ion f or a partially ordered magnetic field. The magnetic 
field is taken to be in the plane of the sky with fractional o rder given by the parameter /„ described in ^3.21 The results of the simulation are plotted as open 
symbols. A fit to the analytical expression given by equation |A3| is plotted as a solid line for each simulation, and the fitted values for e and f appear in the figure 
legend. 

of Strength k = 0.64 was applied and the lower cutoff to the electron energy spectrum was set to 7/ = 19 to generate the observed 
amount of Faraday rotation at each frequency. The observed peak frequency of the component was set to be lypmk = 8.95 GHz 
with a spectral index of a = -0.50. We then ran Stokes / and m; values from this simulation through our analytical peak fitting 
program and found Vpeak = 8.80 GHz with a = -0.49, only 2% less than the input values of the simulation. 
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